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Abstract. We present preliminary results on a search for squarks and gluinos in proton-antiproton 
collisions with a center-of-mass energy of 1.96 TeV and based on 1.1 fb^^ of data collected by 
the CDF detector in the Tevatron Run II. Events with multiple jets of hadrons and large missing 
transverse energy in the final state are studied within the framework of minimal supergravity and 
assuming R-parity conservation. No excess with respect to Standard Model predictions is observed 
and new limits on the gluino and squark masses are extracted. 
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1 Introduction 

Supersymmetry (SUSY) [1] is regarded as one of 
the most promising fundamental theories to de- 
scribe physics at arbitrarily high energies beyond the 
Standard Model (SM). Over the last decades, the Min- 
imal Supergravity (mSUGRA) [5] scenario has been 
extensively studied. In mSUGRA, symmetry breaking 
is achieved via gravitational interactions and the vast 
SUSY parameter space is reduced to only five param- 
eters that determine the low energy phenomenology 
from the scale of Grand Unification(GUT): the com- 
mon scalar mass mo, the common gaugino mass mi/27 
the common soft trilinear SUSY breaking parameters 
^Oi the ratio of the Higgs vacuum expectation values 
at the electroweak scale tan/3, and the sign of the Hig- 
gsino mass term sgn(/i). When R-parityQ|is conserved, 
supersymmetric particles have to be produced in pairs 
and ultimately decay into the lightest supersymmetric 
particle (LSP), usually identified as the lightest neu- 
tralino Xii which constitutes a valid candidate for cold 
dark matter. 

At the Tevatron, the production of squarks (q) and 
gluinos (g), superpartners of quarks and gluons, con- 
stitutes one of the most promising channels because 
of the strong couplings involved. The cascade decay of 
gluinos and squarks into quarks and gluons will result 
in a final state consisting of several jets plus missing 
transverse energy (J^t ) coming from the neutralinos, 
which leave CDF undetected. In this document, we 
present results on a search for squarks and gluinos car- 
ried out on 1.1 fb^^ of data0 collected by the CDF 
detector in Run II. 

^ Email: donofrio@fnal.gov 

^ Rp = (_i)3{B-L)+2s^ where B is the baryon number, 
L is the lepton number, and s is the spin. Rp is 1 for SM 
particles, -1 for supersymmetric particles. 

^ At the time of writing, the search has been updated 
using 1.4 fb~^ of data collected by the CDF detector. 



2 Signal and background samples 

A mSUGRA scenario with Aq ~ 0, sgn(/i)=-l and 
tan/3=5 is assumed. The gluino-squark mass plane is 
scanned via variations of the parameters mg (0-500 
GeV/c2) and mi/2 (50-200 GeV/c^). The pythia [3] 
Monte Carlo program is used to generate samples for 
each mSUGRA point, and ISASUGRA 7.74 is imple- 
mented to predict the SUSY spectrum at the TeV 
scale. Light-flavor squark masses are considered de- 
generate, while 2-to-2 processes involving stop (i) and 
sbottom (6) production are excluded to avoid a strong 
theoretical dependence on the mixing in the third gen- 
eration. Each squark/gluino production subprocess (cor- 
responding to one of the final states gg, qq, qq, qg 
and c.c.) is normalized to next-to-leading order (NLO) 
predictions as estimated using PROSPINO v. 2 [4J. The 
cross sections are computed using CTEQ6.1M PDFs 
and the renormalization/factorization scale is set to 
the average mass of the two final-state supersymmet- 
ric particles. 

SM backgrounds are dominated by QCD multi- 
jet processes where the observed J^t comes from 
partially reconstructed jets in the final state. Monte 
Carlo samples have been generated using pythia 
TUNE A, and data in a low J^t region are used to 
determine the absolute normalization. Other sources 
of backgrounds are Z and W production in association 
with jets, and top and diboson production. For these 
processes, I^t might result from mis-reconstructed jets, 
the presence of neutrinos and muons in the final state, 
or both. ALPGEN v2.1 5 interfaced with parton shower 
from PYTHIA is used to estimate W and Z-|-jets back- 
grounds, taking the normalization from the measured 
inclusive Drell-Yan cross sections [5] . Top and diboson 
processes are estimated using pythia tune a Monte 
Carlo samples normalized to the NLO theoretical pre- 
dictions [7]. Finally, the generated samples are recon- 
tructed using the full CDF detector simulation. 
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3 Event Selection 

Events are required to have a reconstructed primary 
vertex with z-position within 60 cm of the nominal in- 
teraction point, >70 GeV and a tracking activity 
consistent with the energy measured in the calorime- 
ter to reject cosmics and beam-halo background. 
Depending on the relative masses of squarks and gluinos, 
different event topologies are expected. If squarks are 
significantly lighter than gluinos, qq production is en- 
hanced. The squark tends to decay according to q 
qXi, and a dijet topology is favoured, along with I^?t due 
to the two neutralinos in the final state. If gluinos are 
lighter than squarks, gg process dominates. Gluinos 
decay via g qqXiJ leading to topologies contain- 
ing a large number of jets (> 4) and moderate I^^t • 
For Mg « Mq, a topology with at least three jets in 
the final state is expected. Three different analyses are 
carried out in parallel, requiring at least 2, 3 or 4 jets 
in the final state, respectively. In each case, the trans- 
verse energy of the jets must be above 25 GeV, and one 
of the leading jet is required to be central {\r]\ <1.0). 
A minimum azimuthal distance between the jets and 
the I^T is required to reduce the QCD multijets con- 
tribution, and lepton vetoes are applied to reject bo- 
son-|-jets, dibosons and tt backgrounds. 

For each final state, a dedicated study has been 
carried out to define the selection criteria that enhance 
the sensitivity to the mSUGRA signals. The signal sig- 
nificance, S/-\/B, with S denoting the signal and B 
the background number of events, is maximized. Jet 
transverse energies, IjZx and Hrp, defined as the sum 
of the transverse energy of the jets, are the variables 
employed in the optimization. Table 1 summarizes the 
thresholds applied on these variables in the different 
analyses. For the 3-jet analysis, three set of thresh- 
olds, intended for different squark/gluino masses, are 
defined in order to further improve the signal signifi- 
cance across the plane where Mg 
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4 Control regions 

Monte Carlo predictions for SM processes are tested in 
background-dominated regions, referred as control re- 
gions, defined by reversing the selection requirements 
introduced to suppress specific background processes. 
Two different types of control regions are identified: 



Table 1. Set of thresholds employed in the analysis. 
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Fig. 1. T^T and Hx in QCD-dominated samples, after all 
selection cuts for the 4-jet region: data are superimposed 
to QCD®non-QCD background with total systematics un- 
certainties. AH cuts have been applied except the one on 
the variable that is represented. 



— QCD-dominated region: at least one of the leading 
jets is required to be aligned with the I^Zt azimuthal 
direction. 

— EWK/Top-dominated region: lepton vetoes are re- 
versed, and at least one electron or muon, identi- 
fied via the presence of an electromagnetic cluster 
or isolated track, is required in the event. 

Figure [T] shows, for the 4-jet analysis, the IjZt and Hrp 
distributions in the QCD-dominated sample, where 
the shaded band indicates the total systematic uncer- 
tainty on the background estimation (see below) . Good 
agreement is found between data and SM Monte Carlo 
predictions. Similar agreement is observed in all con- 
trol samples considered for each separately analysis. 



5 Systematic uncertainties 

The dominant systematic uncertainty on the signal ef- 
ficiencies and SM background yields comes from the 
3% uncertainty on the jet energy scale. This trans- 
lates into a 10% and 20% uncertainty on the signal 
and background estimations, respectively. Systematic 
uncertainties related to the modeling of the initial and 
final state radiation (ISR/FSR) in mSUGRA signal 
and tt Monte Carlo samples translates into an uncer- 
tainty on the signal efficiency between 3% to 6% de- 
pending on the q/g masses, and in a 8% uncertainty 
on the top background contribution, as defined us- 
ing additional samples generated with modified parton 
showers. Other uncertainties specific for the different 
background processes include: a 10% uncertainty on 
the diboson production cross section due to the un- 
certainties on PDFs and the choice of renormalization 
scale, a 2% systematic uncertainty on boson-|-jets pro- 
duction cross section, and a 10% uncertainty on the tt 
theoretical cross section. Finally, a 6% uncertainty on 
the integrated luminosity is also included. 

Systematic uncertainties on the mSUGRA cross 
sections are dominated by PDFs and renormalization 
scale uncertainties. The Hessian method 8 is used to 
compute the uncertainty due to PDFs. This translates 
into an uncertainty on the theoretical cross section 
that varies between 20% and 34%, depending on the 
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mSUGRA point. The uncertainty due to the renor- 
mahzation and factorization scale is determined by re- 
evaluating the cross section using twice and half the 
nominal scale set in prospino for each production sub- 
process. The resulting systematic uncertainty varies 
between 17% and 23%, depending on the mSUGRA 
point. 



6 Results 

The number of observed and SM expected events cor- 
responding to a total integrated luminosity of 1.1 fb""'^ 
are reported in Table 2, separately for the five analy- 
ses. Good agreement is found between data and Monte 
Carlo predictions. Figure [3 shows the I^?t and Ht dis- 
tributions for the 3-jet(C) and 4-jet analyses. All selec- 
tion criteria have been applied except the one on the 
variable that is represented. Each Figure shows the 
data superimposed on the expected SM background 
and the total systematic uncertainty. For illustration 
purposes, one representative mSUGRA point is shown 
for each analysis region. In case of the 3-jet(C) region, 
a point with Mg « Mq w 340 GeV/c^ is chosen. For 
the 4-jet region, a point with Mg = 287 GeV/c^ and 
Mq = 439 GeV/c^ is considered. 

Since no significant deviation from the SM predic- 
tions is observed, results are translated into exclusion 
limits on gluino and squark production as a function 
of gluino and squark masses. A Bayesian approach at 
95% C.L. is used, where statistical and systematic un- 
certainties are included in the limit calculation tak- 
ing into account correlations between signal and back- 
ground uncertainties. The uncertainties on the NLO 
signal cross sections are also included. 
Figure [3] shows the prospino NLO cross section as a 
funtion of the squark masses for Mg k, 230 GeV/c^, 
and as a function of the gluino mass for Mg k, Mq. 
The yellow band indicates the total systematic uncer- 
tainties on the NLO predictions. 

The figure also shows the curves for expected and 
observed limits, where the crossing point with the nom- 
inal NLO prediction would define the maximum values 
for squark/gluino masses as resulting by this analy- 
sis. The obtained exclusion limits from different squark 
and gluino masses are mapped out in Figure 2] where 
exclusion regions, as determined by other experiments, 
are also presented. 

This search excludes masses up to 385 GeV/c^ at 95% 
C.L. in the region where gluino and squark masses 
are similar, gluino masses up to 280 GeV/c^ for ev- 
ery squark mass, and gluino masses up to 410 GeV/c 
for squark masses below 380 G(N/c\ in a mSUGRA 
scenario with = 0, /x < and tan/3 = 5. 

The results of this analysis in terms of Mg-Mq can 
also be translated into the mSUGRA parameters at 
the GUT scale. Figure [S] shows the excluded regions in 
the mg-m-^y2 plane for tan/3=5, Aq=Q and ^ <0. This 
search improves on the limits from indirect searches as 
determined by LEP2 f9] for values between 75 and 



Table 2. Observed number of data events for the five se- 
lection analyses in 1.1 fb~^, compared with the expected 
events from SM processes. The quoted systematic uncer- 
tainty on the background include 6% uncertainty on the 
luminosity. 
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GeV/c2. 
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Fig. 2. 'EjiT and Ht after all selection cuts for the 3-jet(C) 
(top) and 4-jet (bottom) analyses: data are superimposed 
to QCD and non-QCD background with total systemat- 
ics uncertainties. Distributions for a representative signal 
point is shown for both 3-jet(C) and 4-jet analysis regions. 
In these plots, all cuts have been applied except the one on 
the variable that is represented. 



7 Conclusion 

We have presented preliminary results on searches for 
squarks and gluinos in proton-antiproton collisions with 
a center-of-mass energy of 1.96 TeV at the Tevatron, 
based on 1.1 fb~^ of data collected by the CDF de- 
tector in Run II. In a mSUGRA scenario with Aq = 
0, /i < and tan/3 = 5, we exclude masses up to 
about 385 GeV/c at 95% C.L. in the region where 
gluino and squark masses are similar, gluino masses 
up to 280 GeV/c^ for every squark mass, and gluino 
masses up to 410 GeV/c for squark masses below 



Colliders - SUSY Phenomenology 



Contributed Talk 



CDF Run II Preliminary 



NLO: PROSPINO CTEQ6.1M 
syst. uncert. {PDF © Ren.) 
observed lirr 
expected Mm 



lit 95% C.L. 
lit 95% C.L. 




= 230 GeV/c'' 



CDF Run II Preiiminary 



450 _ 

[GeV/c ] 



L=1.1 fb"* 



NLO: PROSPINO CTEQ6.1M 
syst. uncert. {PDF © Ren.) 
observed limit 95% C.L. 
expected limit 95% C.L. 




g q 



excluded up to 385 GeV/c 



380 GeV/c . In the future, Tevatron is expected to de- 
liver a total integrated luminosity of about 6 fb~^ . This 
will allow either a discovery or further improvement on 
the current limits on squark and gluino masses. 
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Fig. 3. Cross section as a function of gluino/squark masses 
in the case of Mg = 230 GeV/c^ (top) and Mq « Mg (bot- 
tom). The observed and expected limits at 95% C.L. are 
also shown. The yellow band indicates the total systematic 
uncertainties on the NLO prediction. In case of Mg ^ 230 
GeV/c^, squark masses are excluded up to arbitrarily high 



values. In case Mq « Mg, squark and gluino masses are ^— 



Fig. 4. Observed(red region) and expected (dashed-line) 
exclusion limits at the 95 % C.L in the Mq - Mg mass plane. 
Regions excluded by previous experiments are shown, to- 
gether with no-mSUGRA solution region. 
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Fig. 5. Observed(solid line) and expected (dashed line) 
exclusion limits at the 95% C.L. in the (mo,mi/2) plane. 
In the dark-grey region there is no mSUGRA solution. The 
light-grey region indicates the region where m(f ) < ?n(xi). 
The beige and green regions are excluded by LEP2 chargino 
and slepton searches [9], respectively. No mSUGRA solu- 
tion is found in the dark-grey region. The light-grey region 
indicates the region where m(f) < m(x5). The black lines 
are the iso-mass curves for gluinos (horizontal) and squarks 
(elliptic) corresponding to masses of 150, 300, 450 and 600 
GeV/c^ 



